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Abstract 

We investigate gravitino dark matter produced thermally at high temperatures and in decays of a long- 
lived sneutrino. We consider the Non-Universal Higgs Model and a generalized gauge mediation model, and 
in each case identify sneutrino LOSP regions of the parameter space consistent with the mass of the Higgs- 
like boson observed at the LHC. We apply relevant collider and cosmological bounds, including constraints 
from Big Bang Nucleosynthesis and from warm dark matter on large scale structures. Generally, we find 
allowed values of the reheating temperature Tr below lO'"* GeV, i.e. somewhat smaller than the values needed 
for thermal leptogenesis, even with a conservative lower bound of 122 GeV on the Ifiggs mass. Requiring 
mass values closer to 126 GeV implies Tr below 10^ GeV and the gravitino mass less than 10 GeV. 
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1 Introduction 

Of all extensions of the Standard Model of elementary particles, its supersymmetric version (Minimal Supersym- 
metric Standard Model, MSSM) still remains particularly well motivated (for a review, see, e.g., [I]). Several 
mechanisms have been proposed to describe the origin and mediation of the necessary supersymmetry breaking 
in the MSSM, yielding distinctive mass spectra for the supersymmetric partners of the known particles. Among 
these, models equipped with i?-parity predict that the lightest supersymmetric particle (LSP) is stable which 
allows for a possibility that it constitutes dark matter (DM) in the Universe. 

In the MSSM, the most popular DM particles are: the lightest ncutralino, the gravitino (present in the MSSM 
embedded in supergravity) or the axino (in the MSSM extended with a C/(1)pq symmetry) 2 . Different proper- 
ties of these particles require different variants of the history of the early Universe, including the observationally 
determined abundance of DM. The lightest neutralino LSP is considered as perhaps the most natural choice for 
DM as its relic abundance from thermal freeze-out can agree with observations, f2cDM/i^ = 0.113±0.004 3]: for 
a 0(100) GeV, although scenarios in which the lightest neutralinos are produced in decays of heavier particles 
previously dominating the energy density of the Universe, have also been proposed, see, e.g., [4]. 

The abundance of extremely weakly interacting particles (EWIMPs) , like the gravitino or the axino LSP, is 
determined in a way markedly different from that of the lightest neutralino. They are produced in scatterings 
of other particles in the primordial plasma; the abundance of such thermally produced EWIMP is proportional 
to the reheating temperature Tr. In addition to thermal production (TP), EWIMPs can also be produced 
in the decays of the next-to- lightest supersymmetric particles (NLSP), which are usually the lightest ordinary 
supersymmetric particles (LOSP). While in general nearly any MSSM particle can be the LOSP, the most 
natural choices, other than the lightest neutralino, are: the lighter stau (or stop), or the sneutrino. It is this 
last case that we will consider in this paper. 

In non-thermal production (NTP) from LOSP decays, the LOSP lifetime may be long enough for energetic 
decay products to affect the abundances of the light elements (particularly in the case of gravitinos, while much 
less so for axinos). A good agreement of predictions of the standard Big Bang Nucleosynthesis (BBN) with 
observations sets stringent limits on such additional contributions (for reviews, see e.g. [5l[6l[7]). The resulting 
picture can be roughly summarized as follows: the constrains are weakest for light (small energy release in 
decay), short-lived (r < 0(100 sec)) and not too-abundant decaying particles, but above all for decays with a 
small hadronic branching ratio. Therefore, of all possible LOSPs, sneutrino decays are the least constrained by 
BBN since hadronic showers can only be produced in (strongly suppressed) 3- and 4-body sneutrino decays. 
It should also be noted that the present LHC data, albeit quite restrictive for the gluinos and the colored 
scalar partners of quarks of the first and second generation, still allow sleptons and sneutrinos with much 
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lower masses. All this makes supersymmetric models with gravitino LSP and sneutrino LOSP quite attractive 
phenomenologically from the bottom-up perspectiv^. 

In the simplest scenarios of supersymmetry breaking the sneutrino is not the LOSP and there have been just 
a handful of studies devoted to analyzing BBN constraints on sneutrino LOSP [9l [TOl [TTJ [12] in theoretically 
motivated scenarios, such as the Non-Universal Higgs Model (NUHM) [13l [141 [H] or the Generalized Gauge 
Mediation (GGM) model [Ml HH [El [ig . 

In the present letter we re- visit the possibility of gravitino dark matter from both TP and from NTP of 
sneutrino decays in light of recent discovery at the LHC of a Higgs-like boson with a mass of 126 GeV [23 [S]- We 
will show that taking this new result at face value implies a stringent upper bound on the reheating temperature 
Tr below 10^ GeV, and also favors a low gravitino mass region, below 10 GeV. Assuming a conservative lower 
bound on the Higgs boson mass of 122 GeV leads to weaker constraints on and allows a larger gravitino 
mass. 

So far Tr has been allowed to take the largest values for the sneutrino LOSP of all MSSM choices owing to 
the lowest yield at freeze-out. This had a double effect of having the weakest effect on BBN and also requiring 
largest Tr for the TP contribution to make up for the reduced gravitino relic abundance from NTP. However, 
the relatively large Higgs mass implies larger SUSY breaking mass scale, thus also larger masses of gauginos 
whose scatterings dominate TP. This in turn boost the TP contribution (which is proportional to their square) 
and, as a consequence, requires lower Tj^. A larger SUSY breaking scale also implies larger sneutrino mass, 
and, as a consequence, larger yield at freeze-out and more energetic hadronic showers, which translates to more 
stringent BBN bounds. 

In addition to analyzing the effect of the Higgs mass on the scenario, we extend and update previous analyses 
in several directions. First, for small and moderate tan /3 we identify the patterns of soft supersymmetry breaking 
masses at the high scale that can lead to the sneutrino NLSP. In realistic models, sneutrinos are often degenerate 
in mass with right-handed sleptons and/or neutralinos (in addition to degeneracy with corresponding left-handed 
charged sleptons). This leads to many coannihilation channels, which may affect the relic density f^NLSP^^ and, 
as a consequence, also BBN bounds, as they are sensitive to r^NLSP^^- Second, we implement the BBN bounds 
using a state-of-the-art numerical code for solving the relevant Boltzmann equations |22j . As an input parameter 
to that code, we perform a full computation of the hadronic branching fraction of sneutrino decays, including 
quark- antiquark production from on-shell and off-shell electroweak gauge bosons and gauginos. 

The paper is organized as follows. In Section 2, we review the constraints in v LOSP scenarios, presenting 
the assumptions leading to sufficiently light sneutrinos, together with the constraints from the Higgs boson 
mass measurement and BBN. In Section 3, we perform a numerical analysis of v LOSP scenarios, discussing 
the impact of various constraints on the NUHM and the GGM model. We present our conclusions in Section 4. 

2 Review of constraints in sneutrino LOSP scenarios 
2.1 Soft supersymmetry breaking masses at the low scale 

The r-sncutrino, is the lightest of the sneutrinos due to the r- Yukawa coupling driving its mass slightly below 
the sneutrinos of the other two generations, and from now on we will refer to it as simply the sneutrino. The 
sneutrino can become lighter than its charged slepton partner thanks to the electroweak I?-term contributions 
to the slepton and the sneutrino masses. For moderate tan /3 the sneutrino mass after electroweak symmetry 
breaking (EWSB) reads (see e.g. [23]) 

ml^ml + Dl, (1) 

^The case with light sneutrino LSP, constituting even a small portion of DM, is very strongly constrained by direct detection 
experiments [Q; in fact, the results of the XENONIOO experiment require this contribution to be at most of the order of 10~^ of a 
relic sneutrino abundance. 



2 



whereas the mass matrix of the charged sleptons of the third generation is given by 



2 , .-i < Dj^ nir {fi tan 13 - Ar) 



« y mr in tan 13- Ar) m% + Dj^ j' 

In dU) and ([2|) m| and to|; denote soft supersymmetry breaking mass parameters of the superpartners of the 
left-handed and the right-handed leptons, respectively, and Ar stands for the t trilinear parameter, with all 
the parameters evaluated at the EWSB scale. The Z3-term contribution to the sneutrino mass, I?? = ~^AI^, 
is negative, while analogous contributions to the masses of the charged sleptons, Z?| = — ^Af^ and 

D| = M| - M^, are positive. 

The sneutrino mass ((T]) is smaller than the smaller of the eigenvalues of the slepton mass matrix ([2]) if the 
condition 



< - < > — + ^w- T^M'z (3) 

is satisfied. For typical values of the parameters, it follows from the inequality Q that the sneutrino is the 
lightest slepton if m\ < m|, and the left-right mixing in the slepton sector is not too large. For example, with 
n — 1000 GeV and tan/3 ~ 10, the condition ([3]) is satisfied if the splitting between \fm^ and \frn^ is of at 
least about 100 GeV. It is also clear that increasing tan /3 while keeping the other parameters in ([T]) and ([2j 
fixed decreases the mass of the lighter charged slepton of each generation, eventually closing the region of the 
parameters where the sneutrino is the lightest slepton. 

2.2 Conditions for < m|; from renormalization group running from a high scale 

Since sneutrino LOSP disfavors large tan/3, we can, following the method outlined in [24 , obtain solutions of 
renormalization group equations (RGEs) for the soft supersymmetry breaking parameters m|; and m| : 



22" \^ st'''^''^^ 



= "^i,o + CLxMl + CL2Mi + CLQinl^^ + [ i ^ ^ ) + ^2 ^ (5) 



where by to^q (with an additional index 0) for S = E, L,Q,U, D we denote sfermion masses at the high 
scale, while Mi 2 are low-scale U{1) and SU{2) gaugino soft mass parameters. The coefficients cei and CLi 
can be found by solving the 1-loop RGEs, whereas ceih, clq by solving the 2-loop RGEs and identifying the 
leading effects; they are given in Table [T] for some representative choices for the high scale Q and the scale 
ms = ^JWq^fni^ at which electroweak symmetry breaking is evaluatec0. (denoted in literature also as Sq) 
is defined as 

= 5o = tr [FM^^^i^^^ J = m]j^ - m]^^ + tr [m^ „ - 2m2, „ + „ - mi „ + m| q] , (6) 

where m| q are the 3x3 sfermion mass matrices at the high scale, and "rri^H^ are the soft supersymmetry 
breaking masses of the Higgs doublets at the high scale, and gi (51,0) is the U{1)y gauge coupling at the low 
(high) scale. Leading corrections arising due to the r Yukawa couplings are denoted by and 5\ y^; for 

small and moderate values of tan (3 they are small and their only role is to make the third generation of sleptons 

slightly lighter than the first two, but they can become important if the mass parameters ^ / nij^ at the high scale 



^We assume that at the high scales the soft supersymmetry breaking parameters are the same for all three generations; beyond 
that framework, e.g. in models with inverted hierarchy of soft supersymmetry breaking masses, two-loop contributions proportional 

ij^ to values smaller than m|,, 



to squark masses can drive m? to values smaller than mi,, opening up a possibility for yet another example of sneutrino LOSP 



m 
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ms 


CEl CLl CL2 CEU Clq 

Q = IQi" GeV 


CEl CLl CL2 CEU Clq 

Q = 10^6 GeV 


500 GeV 
1000 GeV 


0.47 0.12 0.52 -0.0027 -0.0049 
0.45 0.11 0.51 -0.0026 -0.0048 


0.62 0.15 0.64 -0.0038 -0.0060 
0.59 0.15 0.62 -0.0037 -0.0059 



Table 1: Numerical values of the coefficients cei, cli, cl2, ceu, clq in (jl]) for two representative choices of the 
high scale Q and of the EWSB mass scale ms — ^JfnfTrTi^. 

or the coefficient At in the trilinear coupling of staus are much larger than \/rn^ and \/rn^- In the case when 
the colored particles are much heavier than the sleptons, as is usually the case, one should in principle include 
the leading two-loop contributions to the RGEs in order to obtain 0(10 GeV) accuracy in mass determination. 

Substituting (U) and ^ into ([3]), we see that the sneutrino can be the LOSP in two (mutually not exclusive) 
cases. Note first that in models with = and universal gaugino masses, such as the CMSSM or the NUHMl 
model (i.e. the NUHM with mn^ = mHa)j for which M2 ~ 2Afi, and a high scale > 10^^ GeV, the sneutrino 
cannot be the LOSP, since it is always heavier than the bino. We can then firstly demand < 0, which gives 
m| < Ml, and then the sign difference in the coefficients multiplying in Q and (O can lead to mj^ < m\. 
This possibility is realized in the general NUHM and later we shall also discuss the corresponding mass spectra 
in some detail. The second option is to relax the gaugino mass universality. This possibility is naturally realized 
in GGM models, leading to the sneutrino LOSP and we shall later present some representative examples of 
mass spectra arising in such models, as wellH. 

2.3 Higgs boson mass of 126 GeV 



Recent data from the LHC [201 HI] strongly suggest that the lightest Higgs boson has a mass of approximately 
126 GeV. As we mentioned in the Introduction, this implies a larger supersymmetry breaking scale and implies 
non-trivial consequences for the possibility of having sneutrino NLSP with gravitino LSP. 
At one loop, the lightest Higgs boson mass can be approximated as [35] 



mi w m^z cos^^+ ^-j^ 



m| \ 12to| 



(7) 



where m| is the (defined above) product of the stop masses and Xt — At ~ /i/ tan/3. It is well known (see e.g. 
[571 [2H1 ISni ISI] ) tha-t consistency with the Higgs boson mass measurement at ^ 126 GeV points toward large 
values of ms ^0(1) TeV and values of Xt maximizing the second term in the square bracket in ([7]), with largest 
values achieved for Xt ~ iy/6ms- (The other option of increasing ms so that the logarithmic correction in 
([7]) gives the whole necessary contribution is less natural as it requires very heavy stops.) A solution of the 
one-loop MSSM RGEs gives [21]: 

At = c'^^At^Q - cfj^mi/2 (8) 

« 3c^Ajo-3c^_-^^2Af,omi/2-|-C^/2™l/2 + --- (9) 

^Q « -c^;A2„ + c^^^/2^t■0"^l/2 + cf/2"^?/2 + •■• (10) 



'^c'XAlo + 2c'Xi^2'^t,0mi/2 + Ci/2TO^/2 + 



(11) 



At one loop the values of the numerical coefficients can be expressed as functions of the gauge and top 
Yukawa couplings. In Table [2] we indicate typical values of these coefficients for different choices of Ms and 
gaugino mass patterns. 



^Another way would be to assume large rriQfi, since it would give a negative contribution to However, this would lead 

to large fi, hence would increase the left-right mixing in the stau sector and would thus make the lighter stau lighter than the 
sneutrino. 



4 



coefScient 


At At u u u Q U 
^A ^1/2 ^A '^A,l/2 ^1/2 ^1/2 ^1/2 


range 


-0.4 -2 -0.04 -0.1 -3 -3 1-2 



Table 2: Approximate values of the coefRcients c" in ([8|)- (fTT|) for ms varying from 1 to 5 TeV and two patterns 
of gaugino masses at the high scale Q = 2x 10^^ GeV: universal gaugino case and Mi^o, M2,o <C Afs.o = ™i/2- 



For brevity, in ((9))- (|lip only the terms depending on the high-scale parameters toi/2 and At^ are shown, 
as they suffice for the following argument. From ©-(HI) it is obvious that the easiest way of obtaining a large 
negative Xt is to make the gluino rather heavy; increasing At.o by an equal amount is about five times less 
effective and may threaten to make the stops tachyonic. However, it should be kept in mind that a large mi ^2 
tends to make large; it is of no particular consequence for Xt, as /i enters this quantity multiplied by 1/ tan /?, 
but a large additionally increases left-right sfermion mixing, which, as we discussed in Section [2.11 tends to 
make charged sleptons lighter than sneutrinos (for fixed sfermion masses at the high scale). 

If the soft supersymmetry breaking mass parameters At, roq and mfj are dominated by the RGE contribu- 
tions from the gluino mass, one obtains Xt/ms ~ Xt/ y/rnqmu — —1, which is not very close to the maximal 
stop mixing scenario i2 , optimal for a large Higgs boson mass (the second term in the square bracket in ([7]) 
is — 2 times smaller than its maximal value). A Higgs boson mass of 126 GeV can be then obtained either by 
assuming a rather large 1111/2, or by taking a large negative Atfl, preferably Atfl — —(1 — 3)mi/2 [33l, or else 
by admitting tachyonic stops at high scales [3D], which we shall not pursue further here. 

We are therefore led to the conclusion that a Higgs boson mass of 126 GeV puts an important constraint 
on the possibility of sneutrino LOSP by implying a higher scale of supersymmetry breaking. Lower bounds 
from direct SUSY searches are consistent with this trend but currently not yet as strong. A large Higgs 
boson mass favors large negative At which is usually correlated with large At via RG running; this increases 
left-right mixing in the charged slepton sector and makes a stau lighter than the sneutrino. Such a large At 
most easily originates from a large mi/2 (or a combination of slightly smaller mi/2 and a large negative Atfi), 
which increases fi, thereby also increasing left-right mixing in the charged slepton sector. In the following we 
will illustrate within two SUSY models employing different supersymmetry breaking mechanisms, and both 
allowing sneutrino LOSP, how a large value of mi/2 implied by a heavy Higgs boson leads to strong constraints 
on the reheating temperature resulting from BBN bounds. 



2.4 Bound on the reheating temperature from BBN 

For gravitinos with a mass significantly smaller than the Fermi scale, their present abundance resulting from 
scatterings in thermal plasma [Ml ESI [36] can be approximated by [37] : 

where tuq is the gravitino mass, Mr denote gaugino mass parameters at the low scale and the coefficients 7^. 
can be calculated from 1-loop RGEs for the gaugino masses and gauge couplings: they can be evaluated for 
Tr = 10^ (10^) GeV as 73 = 0.50 (0.67), 72 = 0.51 (0.49), 71 = 0.20 (0.15), for the gluino masses of 900 GeV. 
It is easy to read from this estimate that with mi/2 ^ 1 TeV and m^, = 100 GeV the observed dark matter 
abundance implies a reheating temperature of Tn — 5 x 10^ GeV, which is close to minimal values — 2 x 10^ GeV 
(—2x10* GeV) required by simple models of thermal leptogenesis with zero (thermal) initial abundance of the 
lightest right-handed neutrinos and sneutrinos [38[ 139] . 
The lifetime of sneutrino LOSP can be approximated as 

/ 4 N / \2 /100GeV\"V, "^^ 
TNLSP = 5.9 xloSec —% 1-^^ ' 13 

^ VI GeV/ V "INLSP / \ ^NLSP 
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All scans: mg^ = 500 GcV, toh, ^ 4000 GcV, /i > 



Case 


varied parameters 


fixed parameters 


1 


mo 


mi/2 




-3000 GeV tan/3 = 10 


2 


Ao 


mi/2 


Too 


= 300GeV tan/3 = 10 


3 




tan /3 


Too 


= 300 GeV TOi/2 = 1200 GeV 


4 


'711/2 


tan /3 


Too 


= 300 GeV Ao ^ -3000 GeV 



Table 3: Description of scans over the parameters in the NUHM presented in Figures [T] and [2j 

which can easily be of the order of lO'^ — 10^ sec. For such long lifetimes it is then possible that hadro-dissociation 
processes induced by a subdominant decay process of sneutrino LOSP where a quark- antiquark pair is produced 
can alter the BBN predictions beyond the current observational uncertainties. We shall study this issue in the 
following Section, by numerically analyzing representative examples in two models of supersymmetry breaking 
which allow for sneutrino LOSP. 

3 Numerical analysis 

As we have argued in Section [2?2l models of supersymmetry breaking at the high scale allow a sneutrino LOSP 
only if at least one of the two conditions: = or Mi : M2 : M3 = ai : a2 : as is violated at the high 
scale. (At 1 loop these conditions are renormalization group invariants, hence they can be evaluated at any 
scale.) A violation of the former is manifest in the NUHM while the latter condition can be satisfied in many 
ways. A recently considered scenario is the GGM where it is assumed that the hidden and the messenger sectors 
in models of gauge mediation can be more complicated than what is required in a minimal theory. (The soft 
supersymmetry breaking parameters in both models are given in the Appendix.) 

3.1 The NUHM 

Armed with the above analytical considerations, we will now identify regions of the NUHM parameter space 
where the sneutrino is the LOSP. We will determine if these solutions are consistent with the Higgs mass, 
low-energy observables and early Universe. The details of the scans are given in Table [31 In our numerical work 
we used suspect |40| to solve the renormalization group equations and calculate mass spectra, micrOMEGAs j41) 
for the LOSP relic abundance and Superlso [42] for flavor observables. 

In Figures [Hand [2] we present the LOSP identity and its mass, as well as the mass of the Higgs boson. In both 
panels of Figure [T] and in the right panel of Figure [2] the sneutrino LOSP region is bounded from above at large 
enough values of toi/2. This can be easily understood since to| is a much faster growing function of TO1/2 than 

which is the bino mass squared. At fixed mi ^2 increasing toq, sfermion masses grow and they eventually 
become larger than the bino mass, which explains the bending of the boundary between the sneutrino and bino 
LOSP regions in the left panel of Figure [1] A negative contribution to to|^, which is proportional to D^, has 
to be overcome by some other positive contributions proportional to toq or TOj^/2; otherwise we find unphysical 
regions (marked white) with tachyonic sleptons. It should also be mentioned that the negative contribution to 
mfj (proportional to D^) is larger by a factor of 3/2 than that to to|; the former parameter also receives a much 
larger renormalization group correction proportional to rn'^^^ than the latter. As a consequence, for sufficiently 
small values of toi/2 the lighter stop becomes lighter than the sleptons; the corresponding region of stop LOSP 
is visible in the left panel of Figure [TJ All these effects leads to a lower bound on mi/2; in our scans we find no 
sneutrino LOSP models for < 800 GeV, which, as we shall discuss later, has the important consequences 

for the maximum reheating temperature. As it can be seen in the right panels of Figures [1] and [2J for /i > 
large negative values of ^o result in large off-diagonal entries in the stop mass matrix and lead to very light and 
even tachyonic stops. The appearance of the bino LOSP region in Figure [2] results from the r- Yukawa effect in 
the renormalization group equations: in the leading logarithm approximation, the quantities 5^ ^^^^ in Q and 
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Figure 1: Slices of the NUHM parameter space: mo vs mi/2 (left panel) and Aq vs mi/2 (right panel) with 
the values of mn^ = 500 GeV, — 4000 GeV fixed at the unification scale and tan/3 = 10, ^ > 0. Contours 
of constant LOSP (Higgs boson) masses are shown as dashed (solid) lines. Unphysical regions are marked in 
white. 




Figure 2: Slices of the NUHM parameter space: Aq vs tan/3 (left panel) and mi/2 vs tan/3 (right panel) 
with the values of mo = 300 GeV, m//^ — 500 GeV, niH^ — 4000 GeV fixed at the unification scale and /x > 0. 
Contours of constant LOSP (Higgs boson) masses are shown as dashed (solid) lines. Unphysical regions are 
marked in white. 
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observable 


BR(6 ^ 57) [45, BR(B„ ^ rvr) \M, BR(S^ ^ ^r) [47] AMb, [IS] 


lower bound 
upper bound 


2.8x10-4 0.7x10-'' 0.7x10"^ 12.9 ps"! 
4x10"'' 2.7x10-'' 6.3x10-9 22.5 ps'^ 



Table 4: 95% CL bounds for selected low-energy flavor and electromagnetic observables. Both experimental 
and theoretical errors have been taken into account. 



observable 


D/H 


^He/D 




•^Li/'^Li 


lower bound 


1.2 X 10-^ 




not applied 




upper bounds 


4 X 10-'V5-3 X 10-^ 


1.5 


0.26 


0.1/0.66 (stringent / conservative) 



Table 5: 95% CL BBN bounds based on [35]. The observables are ratios of the element abundances, with the 
obvious exception of Yp which is ^He mass fraction. The lower limit on Yp is irrelevant for constraining the 
abundance of decaying particles. The upper limit on D/H represents a compromise between the commonly used 
average of the best determinations of this quantity and the large spread of the individual results. The use of 
two bounds for ^Li/^Li, a stringent and a conservative one, reflects the uncertainty in estimating the efficiency 
of production/destruction of this element in stellar environment. 

<5i in © can be approximated by 

4... - - -^.vH-^'h. + Al) log (^) , (14) 

where Mgut is the unification scale at which ruHa and yr are evaluated here. As yr increases with growing 
tan/3, we see that moderate values of tan /3 actually help one of the sleptons to become the LOSP. 

Unsurprisingly, for fixed = — 3000GeV, we find a lower bound mi/2 ~ 1 TeV resulting from the require- 
ment that the Higgs boson mass exceeds 122 GeV, the value that we adopt as a conservative lower bound on 
the observable. For fixed mo, the shape of the constant Higgs boson mass contours in the (Ao,mi/2) plane 
agrees with the requirements for maximal stop mixing |33) . For all the points shown in Figures [1] and [2] the 
low-energy observables lie within a conservative 95% CL range quoted in Table HJ As in many unified models, 
the supersymmetric contribution to 6a ^ is too low to explain the observed anomaly [13] . We also checked that 
for all the points of interest squark masses of the first and second generations are well above 1400 GeV, required 
by the LHC data [33]. In the stop LOSP regions, stop masses are often much smaller that ~ 450 GeV which 
is the lower limit from the LHC, but these regions are disfavored anyway, because the Higgs boson mass drops 
below 120 GeV there. 

For the region with sneutrino LOSP shown in the left panel of Figure [T] we calculate the abundances of 
light elements following the method outlined in [22] and apply the observational limits shown in Table [5j A 
representative sample of our results is shown in Figures [3] and |4l We find no constraints for the gravitino masses 
smaller than 7.5 GeV. At tjiq = 10 GeV a part of the parameter space corresponding to mij ^ 500 GeV or, 
equivalently, to to > 10"^ s, is excluded because of too large D/H abundance. For all values of the bounds 
from ^Li/^Li are always more stringent than the D/H bounds. A further increase of rriQ does not change 
this picture much, until the gravitino becomes partly degenerate with the sneutrino, which introduces a strong 
phase-space enhancement of the sneutrino lifetime. This is illustrated in the case with TOq — 250 GeV, for which 
the BBN bounds, while still dominated by ^Li/^Li, become weaker. 

A closer look at the actual predictions for D /H in the considered parameter range of the NUHM reveals that 
even at points consistent with the allowed bounds, the abundance of D is altered with respect to the standard 
BBN value. It is also quite sensitive to the hadronic energy release: if we approximated it as (mp — m^)/3, as 
is often done in the literature, instead of calculating the energy of the qq pair produced in the sneutrino decay, 
then with the conservative ^Li/^Li limit the lower boundary of the respective excluded region in the right panel 
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Figure 3: BBN bounds for the sneutrino LOSP region in the NUHM shown in the left panel of Figure [T] for 
the values of gravitino mass of ttiq = 10 and 20 GeV. For ^Li/^Li the stringent limit was used; the conservative 
limit does not constrain the parameter space. 




Figure 4: BBN bounds for the sneutrino LOSP region in the NUHM shown in the left panel of Figure [T] for the 
values of gravitino mass of = 40 and 250 GeV. For ®Li/^Li the stringent limit was used; the boundary of 
the excluded region with the more conservative constraint for ^Li/^Li is represented by a red dash-dotted line. 
The dotted green line in the left panel shows the change in the lower boundary of the region excluded by D/H 
if a more conservative limit 5.3 x 10^*^ is used [22) . 
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Figure 5: Left panel: BBN constraints shown in the t^; vs nic,Yc, plane for the sneutrino LOSP region shown 
in the left panel of Figured] Dots show the results of our scan with fixed — 2.5, 20 and 250 GeV. Right 
panel: The impact of different estimates of hadronic energy release on the D/H bounds for m^, — 20 GeV. For 
the excluded region marked 'simple Ehad' an approximation E'had — i'm-i' — Wq)/3 was used, while the excluded 
region marked 'full Ehad' corresponds to a computation of i?had involving integration over the full 4-body phase 
space. 



of Figure [3] would shift downwards by as much as ~ 100 GeV (see the right panel of Figure [5]). In other words, 
one would significantly underestimate the sneutrino LOSP region allowed by the constraint. 

In order to understand better the origin of the BBN constraints, we first project all the analyzed points 
onto the Tp vs m^iYp plane; this is shown in the left panel of Figure [5) We also show there the bounds from 
the abundances of those light elements that constrain parameter space regions with sneutrino LOSP. Since 
rip/i^ is roughly proportional to m? (neglecting the opening of additional annihilation channels for increased v 
masses) and since for ttiq <C mp the sneutrino lifetime scales as Tp oc m^m^^, it is easy to understand why, 
with increasing itiq, the constraints from D/H and ^Li/^Li first appear, next tighten up and then eventually 
become weaker. As can be also easily seen from (|13p . a partial degeneracy between ttiq and rrij} causes a much 
larger increase of than the simple power law above implies, hence the BBN bounds become correspondingly 
weaker. Those features can easily be seen for the results of our scan with a three values of fixed gravitino mass 
of ttIq = 2.5, 20, 250 GeV, as a band of dark red dots. 

One may worry that for long sneutrino lifetimes, Tp > 10^ s, the electromagnetic showers produced in 
scatterings of energetic neutrinos from sneutrino decays off neutrinos of cosmic background can affect the BBN 
by altering the '^He/H abundance [9l |49]. In order to verify this we determined that the exclusion plots from 
[5] do not provide additional constraints on our parameter space; we also interpolated the exclusion plots from 
[49] in {mp,Tp, Bh) plane and found no significant impact from '^He/D. 

As we have seen, increasing for a given mp tends to alleviate the BBN constraints. However, for large 
gravitino masses there is another factor that we have to take into account. Non-thermal gravitinos produced in 
sneutrino LOSP decays will have velocities much larger than those characteristic for thermal distribution. Such 
fast moving dark matter particles tend to erase small scales of Large Scale Structures (LSS), especially when 
they constitute a sizable fraction of the dark matter density. Following [501 , we account for these LSS constraints 
by requiring that the root mean square velocity of the non-thermally produced dark matter gravitinos does not 
exceed 1 km/s and that the non-thermal component makes less than 20% of the total dark matter abundance. 
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Figure 6: A comparison of the BBN and LSS bounds for ni^ = 175 and 250 GeV. Long-dashed hnes show 
contours of constant reheating temperature Tr. For ^Li/^'Li the stringent hmit was used; the boundary of the 
excluded region with the more conservative constraint for ^Li/^Li is represented by a red dash-dotted Une. 



The impact of this bound on the NUHM parameter space is shown in Figure |6l where we show the super- 
position of the BBN bound discussed previously and the above LSS bounds for niQ = 175 and 250 GcV. At 
such large ttIq, the LSS bounds become more stringent than the BBN ones (at > 300 GeV the LOSP relic 
abundance, riLOSp'*^: exceeds 20% of the total dark matter abundance, hence iW^^h"^ is also of this order), 
leaving just a small allowed strip in the parameter space. For ttIq > 270 GeV, we find that the LSS bounds 
exclude the entire section of the parameter space that we analyze here. This has important consequences for 
the maximum reheating temperature, since limits on the maximum reheating temperature become weaker with 
increasing gravitino mass. 

A summary of our results is presented in Figure [7] which shows regions in the (mg,,m£;) plane excluded 
by our constraints. It is clear that the BBN bounds alone allow two distinct regions in the parameter space. 
For small itiq < 10 GeV, there are no constraints on mo but the allowed maximum reheating temperature is 
relatively low, T^^^ ~ 10^ GeV. For larger mg, the BBN bounds start constraining the sneutrino mass and 
the maximum reheating temperature increases to ~ 10^ GeV when mc, ~ ttiq. Imposing the LSS bounds closes 
this second region, thus slightly reducing the reheating temperature down to ~ 9 x 10^ GeV. However, now the 
points for which Tr is maximal correspond to Higgs boson masses much smaller than the LHC measurement. 
The requirement that the Higgs boson mass is at least 122 GeV, brings T^'^^ down to 7 x 10^ GeV. 

These bounds on maximum Tr as a function of are shown in the left panel of Figure |8] for the same 
sets of constraints. In the panel we impose the BBN bounds and we show the results with and without the LSS 
bounds and with and without the requirement that the Higgs boson mass is at least 122 GeV. We see that in 
each case the maximum Tr lies close to 10^ GeV, depending on the set of bounds imposed. Without the LSS or 
the Higgs boson mass bounds, this constraint mainly results from the lower bound on as the maximum 

Tr scales roughly as rn~J^. This can be seen in the right panel of Figure [S] where we show the maximum Tr 
versus the Higgs boson mass with and without BBN and LSS constraints. Note that at mh = 126 GeV the 
maximum T^' plunges down as the BBN and the LSS bound become inconsistent with larger values of the Higgs 
boson mass. 

These maximum values of is close to the quoted above lower bound required by simple thermal lepto- 
genesis. It should be noted that the quoted leptogenesis bounds should be treated as indicative rather than 
absolute, since a rather mild mass degeneracy in the right-handed neutrino sector may lower the minimum 
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Figure 7: A summary of the bounds in the NUHM in the {m^,mi}) plane. The thick dashed hne bounds the 
region excluded by BBN, the solid red line marks the boundary of the region excluded by LSS. Thinner dashed 
lines show the maximum reheating temperature, T^^'^, and thinner dotted lines show the Higgs boson mass 
corresponding to Tj^^^. For ^Li/^Li the stringent limit was used; the boundary of the excluded region with the 
more conservative constraint for ^Li/^Li is represented by a red dash-dotted line. 
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Figure 8: Left panel: the maximum reheating temperature as a function of m^, with the BBN, the LSS and 
the Higgs boson mass constraints (m/i > 122 GeV) applied, as well as without one or both of the LSS and 
the Higgs boson mass constraints. Right panel: maximum reheating temperature versus the Higgs boson mass 
without the LSS constraint (upper red dashed line) and with the LSS constraint (lower line). The solid (dotted) 
segments correspond to the cases where the BBN bound is (is not) applied. 



reheating temperature for successful leptogenesis by a factor of a few [STI. 
3.2 GGM models 

Another class of theoretically motivated scenarios in which it is possible to obtain sneutrino LOSP are models 
of Generalized Gauge Mediation. Unlike in the NUHM, where the condition tr(yMg^^j^j.g) = is violated, 
the feature of GGM models that allows for a sneutrino LOSP is a non-universality of the gaugino masses. In 
particular, it follows from ^ that sneutrino LOSP is viable for M2/M1 < 2 at the electroweak scale. We 
shall therefore utilize the freedom of gaugino mass assignment offered by GGM models to reduce M2,o at the 
messenger scale without breaking the universality of the two remaining gauginos, i.e. we shall adopt Mi^o — -^3,0- 
More specifically, we shall assume Mi^ : Af2,o '■ -^3,0 = 5:2:5, which predicts that the lightest gaugino-like 
neutralino is a wino. Eqs. ^ also show that for sneutrino LOSP rnj^ cannot be too large, which for fixed 
gaugino mass scale places upper bounds on the parameters Ai and A2, whose relation to scalar masses at the 
messenger scale is shown in the Appendix. These bounds can be seen in the left panel of Figure [SI which also 
shows that increasing A2 with fixed Ai increases m| q with respect to m|; q, which may lead to a right-handed 
stau LOSP. 

The interplay between the gaugino and scalar mass scales is shown in the right panel of Figure [SJ where 
we keep Ai — 2A2. Increasing gaugino masses while keeping the mass ratios fixed enlarges the range of A2 
(and Ai) for which one of the sleptons is the LOSP. With small values of Ai slepton masses are governed by 
1-loop corrections proportional to gaugino masses and, as follows from (|4]), right-handed staus are the lightest. 
By increasing A one can obtain sneutrino LOSP, because the slepton masses at the messenger scale, q and 

m|;Q, contain contributions proportional to the product of A^ and the square of the respective hypercharge, 
which is larger for right-handed sleptons. 

GGM models do not allow A-terms at the messenger scale. (This can, however, be circumvented by adding 
direct messenger-matter couplings in the superpotential [55], consistently with a Higgs mass of 126 GeV [55].) 
Therefore, in order to have a sufficiently large Higgs boson mass, we have to consider much larger gluino mass 
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Figure 9: Sections of the GGM parameter space: Ai vs A2 (left panel) and Mi vs A2 (right panel) with fixed 
ratio Q : M2fi : Af 3 = 5:2:5 and fixed values of tan/3 — 10, the messenger scale Mmoss — 10^'^ GeV and 
A3 = 20 TeV with /i > 0. Contours of constant LOSP (Higgs boson) masses are shown as dashed (solid) lines. 



M^fi at the messenger scale, which generates radiatively large stop masses and a large negatively At at the 
electroweak scale. The left panel of Figure |9] corresponds to fixed gaugino mass parameters and, therefore, to 
an almost constant Higgs boson mass of 126 GeV; in the right panel, the Higgs boson mass becomes close to 
126 GeV for large values of Mi.o = M^^. At first sight it might seem that one could use the non-universality of 
the gaugino masses to make Mi_o and M2,o much smaller than A/3^0, but a large M^^ also results in a large 
parameter at the electroweak scale, which prevents the sneutrino from being the LOSP. This is also the reason 
for adopting a relatively small value of A3: too large squark masses at the messenger scale also increase /i. All 
this results in sneutrino LOSP masses of about 1 TeV, which is much larger than in the case of the NUHM. 

Models with gauge mediation of supersymmetry breaking have the advantage that the leading contributions 
to the soft masses are flavor-diagonal, while the subdominant gravity- mediated contributions, of the order of 
tHq, do not have to exhibit any such structure. This leads to a natural suppression of the FCNC's, but also 
has the obvious consequence that TOq ^ mo , with the precise hierarchy depending on details of an appropriate 
flavor model. For this reason we do not consider larger than 20 GeV, hence the reheating temperature 
consistent with the measured dark matter abundance is much lower for GGM than for the NUHM. Nonetheless, 
we find some BBN constraints for m^, ~ 0(10) GeV; they are shown in Figure [TOl 

4 Conclusions 

In this paper, we have analyzed models of gravitino dark matter with the r-sneutrino as the lightest ordinary 
supersymmetric particle. We have shown that if the scale at which supersymmetry is broken is close to the 
unification scale, the sneutrino can be the LOSP either if gauginos are non-universal or < at the high 
scale. We have then performed a detailed study of representative examples of these two possibilities: one arising 
in the NUHM and the other in models of generalized gauge mediation of supersymmetry breaking. 

We have calculated the changes in the BBN predictions due to hadronic showers from sneutrino decays 
D vGqq, calculating the hadronic energy release by a numerical integration over the phase space of the 
produced particles. In the cases in which the D/H bound provided the strongest constraint, we observed big 
changes of the excluded regions of the parameter space between this calculation and one using a simplified 
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Figure 10: BBN bounds for the sneutrino LOSP region in GGM model shown in the right panel of Figure [5] 
for values of gravitino mass = 10 and 20GeV. For ^Li/'^Li the stringent limit was used; the boundary of 
the excluded region with the more conservative constraint for ^Li/^Li is represented by a red dash-dotted line. 



formula for the hadronic energy release. 

We showed that in models of gravitino dark matter with sneutrino LOSP, the four classes of constraints that 
can be applied: (i) the BBN constrains, (ii) the constraints on the large structure formation due to a presence 
of free-streaming decay products of !>, (iii) the Higgs boson mass bounds derived from the LHC data and (iv) 
the bounds on reheating temperature required by simple models of thermal leptogenesis are inconsistent, albeit 
the maximum reheating temperature is only 2 — 3 times smaller than the value suggested by the leptogenesis 
bound. (This was clearly visible in the NUHM; with gauge mediation, such large reheating temperatures were 
unattainable due to theoretical constraints on the gravitino mass.) Therefore, our results challenge the notion 
that models of gravitino dark matter with sneutrino LOSP are compatible with simple thermal leptogenesis. 
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Appendix 

Here we collect expressions for the soft supersymmetry breaking parameters at the high scale in the NUHM 
and the GGM model. 

The parametrization of the NUHM is very simple 

Mifi = M2^o = ^^3,0 = ™i/2 (15) 
2 2 2 2 2 2 

while TTijj^, iTT-Ha^ ^0, tan/? and sgn(/i) can are free parameters. 
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In GGM models, the soft supersymmetry breaking masses at the high scale in the notation of [19] read 

Mrfi = (aj47r)Ai for r = 1,2, 3 (17) 

"^Q.O = (8/3)(a2/16^2)^2^(3/2)(«2/^g^2)^2^(^/30)(^2/^g^2)^2 (18) 

ml, = {8/3){al/167T^)Al + {8/l5){al/16n')Al (19) 

ml^, = {8/3){al/167T^)Al + {2/l5){al/16n^)Al (20) 

<.o = m%^^^ml^ = {3/2){al/16n^)Al + {3/W){al/m7r')Al (21) 

<,o = (6/5)(a?/16^^)A2. (22) 

The trilinear scalar couplings are all equal to zero and tan /3, sgn(/i) are free parameters. 
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